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ABSTRACT 

We present the temperature power spectra of the cosmic microwave background (CMB) derived 
from the three seasons of data from the Atacama Cosmology Telescope (ACT) at 148 GHz and 218 
GHz, as well as the cross-frequency spectrum between the two channels. We detect and correct 
for contamination due to the Galactic cirrus in our equatorial maps. We present the results of 
a number of tests for possible systematic error and conclude that any effects are not significant 
compared to the statistical errors we quote. Where they overlap, we cross-correlate the ACT and the 
South Pole Telescope (SPT) maps and show they are consistent. The measurements of higher-order 
peaks in the CMB power spectrum provide an additional test of the ACDM cosmological model, 
and help constrain extensions beyond the standard model. The small angular scale power spectrum 
also provides constraining power on the Sunyaev-Zel'dovich effects and extragalactic foregrounds. We 
also present a measurement of the CMB gravitational Icnsing convergence power spectrum at 4.6(t 
detection significance. 

Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

The current generation of arcminute resolution cosmic 
microwave background (CMB) experiments is provid- 
ing researchers with a precise view of CMB anisotropics 
over a range of scales (500 < £ < 10000). Over 
the so-called Silk damping tail of the CMB (500 < 
£ < 3000) these observations are revealing the sub- 
tle effects that inflationary physics, primordial he- 
lium density and the energy density in relativistic de- 
grees of freedom have on the acoustic oscillations in 
the photon-baryon plasma in the radiation-dominated 
era. Rapid progress in measurements of the damping 
tail of the power spectrum has been achieved over a 
span of a few years by a number of exp eriments, no- 
tably the Cosmic Background Imager (CBIi lSievers et all 
|2009( ). Arcminute Cosmology B o lomet er Array Re- 
ceiver (ACB AR: iReiehardt et"ni I2009D. QUEST a t 
DASI (QUaD; iBrown et al.ll2009HFriedman et al.ll2009D 
the Atacama Cos mology Telescope (AC T; iSwetz et al.l 
2010HFowler et al.ll2010l : lDas et al.|[2011bl hereafter Dll: 
Dunklev et al.ll201lD and the South Po le Telescope (SPT; 
Keisler et all l20lH iStorv et al.l |2012[ ). Minute distor- 
tions of the CMB anisotropics due to the gravita- 
tional lensing of CMB photons by large-scale structure 
have now been detected using CMB data alone, both 
in percent-level al teration of the damping tail acous- 
tic pe ak structure (jPlll . iKeisler et al.l 120111 IStorv et al.l 
I2012t ). and in the subtle no n- Gaussian signa ture it in- 
duces i n the statistics of CMB anisotropics ()Das et al.l 
I2011at Ivan Engelen et al.l I2012D . When added to the 
cosmic variance-limited CMB power spectr um measure- 
ment s ate< 1000 by the WMAP sateUitc (|Larson et al.l 
120 111 ), the damping tail measurements are providing ad- 
ditional dynamic range resulting in improved constraints 
on inflationary parameters such as the tilt and run- 
ning of the primordial power spectrum. On smaller 
scales {i > 3000) the primordial CMB signal dimin- 
ishes and emission from radio galaxies and dusty star 
forming galaxies, as well as the thermal and kinematic 
Sunyaev-Zel'dovich effects (jSunvaev fc ZeldovichI I1972D 
arising from the scattering of CMB photons by hot gas 
in galaxy clusters, dominate the power spectrum. Along 
with the damping tail measurements, rapid progress has 
been made in the measurements and modeling of this 
high-multipole tail (iLuekcr et aLI 120101: IShirokoff et al.l 
[201lt [dTTF iDunklev et al.ll201lL iReichardt et al.ll2012i r 
These measurements have been used to estimate the 
thermal and kinematic SZ contributions to the power 
spectrum, as well as to model the cosmic infrared back- 
ground (CIB) power spectrum arising from th e mm- wave 
emiss io n from dusty h igh re d shift galaxies (iHall et al.l 
20T0t [P^klcv et al.l I20TTI: lAddisqn et al .1 I2012bllal: 
Haiian et al.l l2012t IReichardt et all 120121 : IZahn et al l 
20121) . ^ 

In this work, we present the measurement of the power 
spectra of CMB temperature anisotropics at 148 GHz 
and 218 GHz from a subset of ACT observations per- 
formed over the 2008, 2009 and 2010 observing seasons, 
and covering approximately 600 deg'^ of the sky. This is 
approximately twice the survey area used in a previous 
measurement of the ACT power spectrum ljDllD . Addi- 
tionally, we also present an updated measurement of the 
gravitational lensing power spectrum from the equatorial 



TABLE 1 

Observations used in power spectrum estimation 



Year Key R.A. Range Dec. Range Area rip^ 

(deg2) 



Equatorial Strip ACT-E 

2009 3e 20''40™ .. 3'^20™ -1?5 .. 1?5 300 5 

2010 4c 20''40™ .. 3'^20™ -1?5 .. 1?5 300 5 

Southern Strip ACT-S 

2008f 2sf 00'^22™ .. OQ^HT^ -55?0 .. -50?0 292 4 

2008 2s 041^08™ .. 07'^08'" -55?2 .. -51?2 146 2 

2009 3s 041^08™ .. 07'^08™ -55?2 .. -51?2 146 2 

2010 4s 04^08™ .. 07'^08™ -55?2 .. -51?2 146 2 

number of patches 

strip. iDunklev et al.l (|2013D use the temperature band- 
powers reported in this work to generate likelihood func- 
tions which form the basis of the cosmological parameter 
constraints reported by Sievers et al. (2012). 

The paper is organized as follows. In Section |2| we 
describe the observations, the generation of maps from 
time ordered data, and the estimation of the beam trans- 
fer functions. We discuss the calibration of the maps in 
Section |3| Section S) describes the pipeline used to pro- 
cess the maps into the angular power spectrum. Treat- 
ment of point sources and other foregrounds is discussed 
in Section |S] Simulations used to test and validate var- 
ious portions of the power spectrum estimation pipeline 
are described in Section |6| The power spectrum results 
and consistency checks are presented in Section [71 and 
the CMB lensing results are discussed in Section [5] We 
conclude in Sectional 

2. OBSERVATIONS AND FIELDS 

ACT is a 6-meter off-axis Gregorian telescope situ- 
ated in the Atacama desert in Chile at an elevation 
of 5190 m. act's Millimeter Bolometric Array Cam- 
era (MBAC) has three channels operating at 148 GHz, 
218 GHz and 277 GHz . The instrument is described in 
detail in ISwetz et"aI1 (|201lD . Between 2007 and 2010 
ACT observed mainly along two constant-declination 
strips on the sky: one running along the celestial equator 
(hereafter the equatorial strip or ACT-E), and the other 
along declination -55° in the southern sky (the south- 
ern strip or ACT-S). The observations were performed 
over Nov 8 ~ Dec 15, 2007; Aug 11 - Dec 24, 2008; May 
18 - Dec 18, 2009, and Apr 6 - Dec 27, 2010. Here 
we present the angular power spectrum measurements 
of 148 GHz and 218 GHz observations from 300 deg^ on 
ACT-E and 292 deg^ on ACT-S. Table [1] summarizes the 
observations from various seasons that were used in esti- 
mating the power spectrum, and also defines shorthand 
notations for maps, e.g., the 2009 season equatorial map 
is referred to as the "season 3e" or simply the "3e" map. 

2.1. Equatorial Observations 

Observations on the ACT-E strip were performed in 
the 2009 and 2010 seasons, and run along the celes- 
tial equator with a right ascension span of 100 degrees, 
and a width of 3 degrees along the declination direc- 
tion. For the power spectrum analysis, we make sin- 
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Fig. 1. — Equatorial maps (ACT-E) made from 2009 (upper panel) and 2010 (lower panel) 148 GHz observations filtered to emphasize 
modes in the range I = 500 — 2500. The four data splits in either season were co-added to make this plot. Also delineated are the patches 
used for computing power spectra. 



gle season maps, and following iDllI we divide the data 
within each season into four splits in time, by distribut- 
ing data from roughly every fourth night into a different 
split, generating four split-maps, each of which is prop- 
erly cross-linked. The maps are also spatially divided 
into five patches on which power spectrum estimation is 
performed separately. We explicitly avoid the edges of 
the maps where the cross-linking is poor and the noise 
is inhomogeneous. A representation of the season 3e and 
season 4e 148 GHz maps and patches are shown in Fig.[TJ 
The two seasons share the same footprint on the sky, 
and common patches were defined to facilitate the com- 
putation of cross season power spectra. Fig. [2] shows 
the noise power spectra of the ACT-E maps by season 
against the CMB-only theory. For most seasons, and 
for 148 GHz, on largest angular scales (£ < 500) atmo- 
spheric noise dominates, while for intermediate angular 
scales (500 < t < 2500) fluctuations in the CMB dom- 
inate the variance. At smaller angular scales detector 
noise becomes the most significant contribution. 

2.2. Southern Observations 

The observations made on the southern sky across var- 
ious seasons had different footprints, requiring a some- 
what involved strategy for efficiently computing the 
power spectrum. Filtered versions of various season maps 
are shown in Fig. [31 The largest coverag e was obtained in 
the 2008 season (the same area on whichlmil was based) . 
When computing the power spectrum within the 2008 
data set, we used four large patches collectively cover- 
ing 292 deg^ (we refer to this full footprint as "season 
2 south full" or season 2sf in short). For computing the 
power spectrum within the other two seasons, as well as 
to compute the cross-power spectra between any pair of 
the three seasons it was necessary to define another set 
of two patches (shown by the smaller contiguous rectan- 
gles in Fig. [3]) that had a common footprint across the 
seasons. As discussed in Section |4l care was taken not to 
double count information while combining the different 
spectra. The noise power spectra of each of the season 
maps are displayed in Fig. U against the CMB-only the- 
ory. The season 3s map is mostly noise dominated on 
all scales in either frequency - we keep this season in 
our analysis to tease out information from cross-season 
spectra, but the season 3s-only spectrum is heavily down- 
weighted in our likelihood. 

2.3. From Time- Ordered Data to Maps 
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Fig. 2. — Noise spectra for each season for the ACT-E maps 
for 148 GHz (upper panel) and 218 GHz (lower panel). The red 
solid line shows the CMB-only spectrum. At 148 GHz the power 
spectrum is sample variance limited at £ < 2500, while at 218 GHz 
detector and atmospheric noise dominate on most scales. 



Details of the map-making pro cedure starting from 
time-ordered data can be found in lDiinner et all (|2012D . 
T he maps were produc ed using the algorithm described 
in lDiinner et all ()2012D . and the 148 GHz 2008 data are 
identical to the maps presented there. We present a 
short summary of the mapping procedure here. Before 
mapping, we reject any detector timestreams that have 
too many spikes (such as cosmic ray hits) identified in 
the data. The cuts were calculated with differing sen- 
sitivites in the spike finder for different bands and sea- 
sons; the threshold for each band/season is set to corre- 
spond to roughly the same fraction of data rejected be- 
cause of spikes. The thresholds are 11/9/6 spikes per 10- 
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Fig. 3. — Southern maps (ACT-S) made from 2008 (top panel) 2009 (middle panel) and 2010 (bottom panel) 148 GHz observations 
filtered to emphasize modes in the range I = 500 — 2500. The four data splits were co-added to make this plot. Also delineated are the 
patches used for computing power spectra. The smaller two patches common between the three maps are used to compute cross-season 
cross-power spectra. The four larger patches for season 2sf are used to compute the full footprint 2008-only cross-power spectrum. Areas 
of large noise or stripes arc heavily down weighted in the analysis. The color scale is the same as Fig[T] 
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Fig. 4. — Noise spectra for each seasons for the ACT-S maps for 
148 GHz (upper panel) an 218 GHz (lower panel). The red solid line 
shows the CMB-only spectrum. Season 3s is significantly noisier 
than the other two seasons. Note that the combination of se asons 
3s a nd 4s is more sensitiv e than season 2s which was used in ID 111 
and lDunklev et al.l 120111) . 
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94 GHz map l|Jarosik et al.|[2011l ) (top) with the map of the same 
region made from ACT 148 GHz (bottom) observations (co-added 
across seasons). All maps have been high-pass filtered with a cos^ £- 
like filter that goes from to 1 for 100 < I < 300. Agreement 
between the CMB features in the two maps is clear by eye. 

15 minute timestream for the 218 GHz 2008/2009/2010 
data, and 9/11/11 for the 2008/2009/2010 148 GHz data. 
We then interpolate across gaps in the remaining detec- 
tor timestreams and deconvolve the effects of the detec- 
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Fig. 6. — Side by side comparison between the ACT map (co-added across seasons) and the SPT map for the same region of the sky. 
The left panel shows the ACT map high-pass filtered with a cos^ £-like filter that goes from to 1 for 100 < £ < 300, and the center 
and r ight panels show the ACT and SPT maps respectively under the same high-pass filter used in the SPT data release llSchaffcr ct al.l 
120111 ). Agreement between the CMB features in the two maps is clear by eye. It is noteworthy that the instrumentation, scan strategy, 
and analysis methods for these two experiments are completely different. 



tor time constants and a (known) filter from the readout 
electronics. Next we remove an offset from each detector 
and a single slope common acros s the array. We the n 
estimate the noise as described in lDiinner et al.l (|2012| ). 
using a model that finds correlations across the array, 
and measures the power spectra of those correlations in 
frequency bins and the power spectra of the individual 
detectors after the correlations are removed (the domi- 
nant correlated signal is a common-mode atmosphere sig- 
nal, but both higher order atmosphere signals and elec- 
tronic noise produce correlated noise across the array). 
We then use a preconditioned conjugate-gradient (PCG) 
iterative algorithm to solve for the maximum-likelihood 
maps. At the same time we solve for the values of the 
timestreams in regions where the data have been cut out. 
The cut samples are assumed to be decoupled from the 
sky; the solution is effectively using the noise model to 
interpolate across the gaps. We do this procedure twice; 
the second time we subtract off the first solution from 
the timestreams to avoid any signal in the data from bi- 
asing the noise estimation. We find that the maps are 
typically unbiased to better than one part in 10"'^ and in 
all cases the transfer function is much smaller than the 
beam error in the i ranges we use for science and calibra- 
tion. Simulations show that the maps typically converge 
within a few hundred PCG iterations. 

2.4. Beam Transfer Functions 

The beam s are estimated indepen dently for each array 
and season (jHasselfield et al.]|2013l in preparation) from 
observations of Satu rn following a proce dure similar to 
the one described in iHincks et al.l (|2010( ) . Radial beam 
profiles from the planet maps are transformed to Fourier 
space by fitting a set of basis functions whose analytic 
transform is known. The fit yields the beam transform as 
well as a covariance matr ix following a procedure similar 
to that discussed in iDlll . The transform is subsequently 
corrected for the mapper transfer function, the solid an- 



gle of Saturn, and the difference in Saturn's spectrum 
compared to the CMB blackbody spectrum. Because any 
location in the ACT CMB maps contains data from many 
different nights, the effective beam in the maps is broad- 
ened relative to the planet-based beam due to pointing 
variation from night to night. This pointing variation 
(~ 6") is modelled as having a two-dimensional Gaussian 
distribution, and the standard deviation is measured by 
comparing the shape of the beam obtained from stacked 
radio sources to the planet-based beam transform. The 
error in the beam due to the pointing correction is in- 
cluded in the final beam covariance matrix. The covari- 
ant error in the beam is obtained after fixing the normal- 
ization of the beams at £ = 700 (1500) for the 148 GHz 
(218 GHz) array. The calibration error is thus separated 
from the covariant error due to beam shape uncertainty, 
which is by construction a.t £ ~ 700 (1500). 

2.5. Comparison with WMAP and South Pole Telescope 

Maps 

Fig. [5] shows a comparison betwee n the WMAP 7- 
year 94 GHz map (jJarosik ct al.ll2011| ). and the 148 GHz 
and 218 GHz ACT maps on the same region of the sky. 
On the one hand, it exemplifies how our map-making 
pipeline faithfully reproduces all the large-scale CMB 
features seen in the WMAP map, and on the other hand 
it portrays the significantly higher resolution afforded by 
ACT over WMAP. This figure is a visual representation 
of the fact that the transfer function in our maps are 
unity down to small angular scales {£ ~ 300) or large 
spatial scales (0?6) - this proves highly beneficial for cal- 
ibrating our maps against the WMAP maps, as discussed 
later. 

We also compare ACT maps with pu blicly released 
maps from the South Pole Telescope (jSchaffcr et al.l 
\20m in the region of overlap in the southern strip. Fig.[B] 
shows a side by side comparison between the ACT map 
co-added across all seasons and the SPT map on the same 
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Fig. 7. — Top: Two-dimensional cross-correlation function in 
real space between WMAP and ACT 148 GHz map (left panel), 
between ACT 148 GHz and ACT 218 GHz maps (middle panel), 
and between the AC T 148 GHz and the SPT 150 GHz map from 
ISchaffer ct"an iWui) . Note that the ACT-SPT cross-correlation is 
anisotropic as the SPT map is filtered in the direction to sup- 
press modes below £x < 1200. Bottom: One-dimensional binned 
version of real space correlation functions. The blue dots represent 
the ID ACT 148 GHz X WMAP correlation, while the green trian- 
gles represent the ID auto-correlation of ACT 148 GHz maps after 
the SPT filter has been applied to them. Finally, the red squares 
denote the ID cross-correlation of the ACT 148 GHz map with the 
SPT map. 

region of the sky. In the middle panel, the ACT map has 
been filtered with the same filter used for the SPT map. 
The similarity between the two maps is clear by eye, and 
speaks to the high quality of both m easurements. A com - 
parison of this figure with Fig. 13 of lDiinner et al.l (j2012D 
(whieh used only 2008 observations) shows the improve- 
ment in the noise properties of the ACT map from the 

combination of multi-season observ ations. 

Using the techniques described in lHaiian et al.l ()2011h . 
we also studied the two-dimens ional real-spaee eor rela- 
tion function (cf. equation 3 of iHaiian "eFalllMl be- 
tween the ACT and the SPT maps. Both the 2D cross 
correlation and the binned ID version are shown in Fig. [71 
The cross-correlation is anisotropic as the SPT map is 
filtered in the £x direction to suppress modes below 
£x < 1200, but the overall agreement between the two 
experiments is excellent. 

3. CALIBRATION 

The final map calibration is performed in two stages: 
first, the 148 GHz map from the lowest-noise season is 
calibrated against the WMAP sky map, and then the 
148 GHz maps from the remaining seasons and 218 GHz 
maps from all seasons are calibrated against the WMAP- 



calibrated 148 GHz map. 

3.1. WMAP Calibration 

ACT map-making and observing strategies result in 
maps with unbiased large-angle modes that can be com- 
pared to WMAP maps of the same region. The maps 
are cross-linked, i.e. every point in the survey has 
been observed during both its rising and setting. The 
cross-linke d data are fed to a m ap-making pipeline de- 
scribed in iDiinner et al.l (j2012f ) that allows the recon- 
struction of all modes in the map without biasing the 
large-angle modes. The transfer function of the maps 
is unity to better tha n 1% at angular s cales corre- 
sponding to £ > 300 piinner ct al.ll2012D . We cali- 
brate the 148 GHz ACT maps d i rectly to WMAP 7- 
year 94 GHz maps ([Jarosik et al.l 120 111 ) of the identi- 
cal regfonsjrsmgthecross-correlation method described 
in IHaiian et"aI1 (|201lD . By matching the ACT- WMAP 
cross-spectrum to the ACT po wer spectrum and the 
WMAP 7-year power spectrum (jLarson m 
the range 300 < £ < 1100, we calibrate the 148 GHz 
ACT spectr um to 2% fractiona l temperature anisotropy 
uncertainty (jHaiian et al.ll201lD . Calibration to WMAP 
is done on the deepest seasons for both ACT-E and ACT- 
S strips, which correspond to season 4e (2010 observing 
season) and season 2sf (2008 observing season) data re- 
spectively. These calibrated maps are used as references 
to calibrate other seasons as described below. 

3.2. Relative Season Calibration 

Once the deepest 148 GHz season has been calibrated 
with respect to WMAP, we cross-correlate that map with 
a 148 GHz map from another season, take the ratio of the 
cross-season cross-power spectrum to the in-season cross- 
power spectrum, and fit for a calibration factor. For ex- 
ample, on the equator, the 148 GHz season 4e map is 
calibrated with respect to WMAP. We then compute the 
ratio (74ex4e^^3ex4e estimate the calibration for the 
season 3e 148 GHz map. This internal method lets us 
use a much wider range of angular scales (we use an £ 
range of width 2000 starting at £ = 500 for 148 GHz and 
£ = 1000 for 218 GHz) than possible with WMAP. Using 
this method, we achieve the following relative calibra- 
tion uncertainties (expressed as ax-Y for season X cali- 
brated against season Y) : cr^e-Ae — 0.7%, a3s-2s — 3%, 
and CT4s_2s = 3% for 148 GHz, and a^e-Ae — 2%, 
o'3s-2s — 9%, and a4s-2s = 4% for 218 GHz maps. Note 
that for the ACT-S season 3s maps the calibration un- 
certainties are higher, as is expected from the fact that 
this season is largely noise dominated on most scales (see 
Fig. [5]). The internal spectrum from this season gets 
highly downweighted in our likelihood. To tie together 
the 148 GHz and 218 GHz internal calibrations, we finally 
calibrate the best season 218 GHz map with respect to its 
148 GHz counterpart, achieving cr2s(2i8)-2s(i48) — 1-3%, 
and cr4e(2i8)-4e(i48) — 1-7%. This gives the overall cali- 
bration of the reference 218 GHz map as 2.4% for ACT-S 
and 2.6% for ACT-E. 

4. TEMPERATURE POWER SPECTRUM ANALYSIS 

The power spectrum a nalysi s methods used here are 
essentially the same as in iDlll . the major modifications 
being due to the smaller extent of the ACT-E maps in 
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the declination direction compared to ACT-S, and the 
muhi-season nature of the spectra. 

4.1. Preprocessing of maps 

We foUow lDllI and apply a high-pass filter to the maps 
to suppress the large angular scale modes {£ < 500) which 
are not as well constrained as others, and can bleed power 
into the smaller-scale modes. Next, we prew hiten the 
maps through real space operations as in ID 111 . 

4.2. Data Window 

Each patch is then multiplied with a data window be- 
fore the power spectrum is computed. The window is a 
product of three components: a point source mask, an 
apodization window, and the riobs map giving the num- 
ber of observations in each pixel. The point source mask 
is further described in Section [521 To simplify the appli- 
cation, we create a single nobs map per patch by adding 
the individual nobs maps from all the splits involved (four 
splits for the single frequency spectrum, and the 8 splits 
for the cross- frequency or cross-season spectrum), and 
apply this as a weight function. This essentially down- 
weights the poorly observed regions of the patch. An 
additional step is applied to the ACT-E patches. Since 
the ACT-E strip is only 3 degrees wide, the absolute 
Fourier space resolution in the declination direction is 
A£y = 120. This leads to instability in the mode-mode 
coupling calculations due to poor sampling of power in 
the Fourier space. To remedy this, we extend the patches 
in the declination direction by adding a 0.7-degree- wide- 
strip of zero valued pixels on either side, such that the 
final declination width of the zero-padded patch is 4.4 de- 
grees. To minimize ringing from the edges of the patch, 
we also apply an apodization window which is generated 
by taking a top-hat function that is unity in the center 
and zero over 10 pixels at the edges of the original patch, 
and convolving it with a Gaussian of full width at half 
maximum of 2.'5 for ACT-S and 14.'0 for ACT-E. An- 
other addition to the pipeline for the ACT-E patches is 
the application of a Galactic dust mask (see Section [5T|) . 
Monte Carlo simulations demonstrate that we retrieve an 
unbiased estimate of the spectrum with these additions 
to our well-tested pipeline. 

4.3. Binning of the power spectrum 

It is important to note that the narrow inherent width 
of the ACT-E strip, as well as the smaller dimensions of 
the mulitseason ACT-S patches prompted us to adopt 
wider bins for the bandpowers than were used in ID 111 . 
Most notably, over the acoustic peaks {£ < 2000) the 
bi ns us ed have a width A£ = 100 as opposed to A£ — 50 
of iDlll . This choice is motivated by the fact that with 
the finer binning, adjacent bins remain significantly cor- 
related for the ACT-E spectrum. In addition, evaluating 
the covariance using full end-to-end Monte Carlo simula- 
tions is prohibitively expensive given the iterative nature 
of our map-making process. Conversely, more tractable 
approximations might not be good enough to provide the 
precision deserved by the high quality of the data. With 
the larger bins we have verified that the bin-to-bin corre- 
lations never exceed 10% and are much smaller than 10% 
for most bin pairs, allowing us to treat the bandpowers 
as statistically independent (cf. Section m|) . 



4.4. Cross-Season Cross- Spectrum Estimation 

To obtain unbiased estimates of the cross spectrum 
we follo w the same steps as enumerated in Section 3.6 
of iDlll . which involve deconvolving a mode-coupling 
matrix that accounts for the effects of beam profile, 
prewhitening, filtering, pixelization, and windowing. 
For same-season cros s spe ctra, the combinatorics are 
exactly the same as in iDllI : we compute 6 cross spectra 
per patch for the single-frequency spectrum (the 4 
"splits" in time giving rise to the 6 cross spectra), 
and 12 for the cross- frequency spectrum (avoiding 
crossing the same splits that contain data from the 
same nights). For cross-season spectra we combine 
all 16 cross-season cross-spectra, as each split from 
one season has independent noise from the splits in 
the other season. For each frequency and season 
pair, the cross-spectra from the patches are combined 
with inverse variance weighting. This results in a 

set of three bandpowers {cf >^3e^ (^3ex4e^ ^4ex4e| 

for ACT-E and a set of six bandpowers 

2sx2s ,^2sx3s /^2sx4s /o3sx3s /-^3sx4s ^4sx4s\ r 
b ^^b ^^b ''-'6 ^^b y^b ) '^"l 

ACT-S, for each of the two same-frequency pairs 
148 X 148 and 218 x 218. For the cross-frequency 

spectra 148 x 218, where cj^'^'^^^-'^^'^*-^^^'' is distinct from 
^4e(i48)x3e(2i8)^ ^ ^^^^ bandpowers for 

ACT-E and nine bandpowers for ACT-S. These add up 
to a total of 10 cross-power spectra for ACT-E and 21 
cross-power spectra for ACT-S that enter the likelihood 
separately with their individual bandpower covariance 
matrices. Additionally, there are six cross-power spectra 
coming from the full- footprint 2008 ACT-S map (2sf), 
which is added, with proper attention to the overlap 
between the s2f and s2 patches. 

4.5. Bandpower Covariance 

For each cross-power spectrum above, we evaluate a 
bandpower covariance matrix: 

Q{aAxi3B);{-tCXTD} _ ^^^aAxpB _ ^^oAx^JS^^j 

x(c,-^x^--(cr^^-)))(l) 

where we used Greek indices a, /3, etc. to denote the 
seasons 3e, 2s, etc. and the uppercase Roman numer- 
als A, B, etc. to denote frequencies. The analytic ex- 
pression of the general term of this covariance matrix 
Q{aAxi3B).{jCxTD) discussed in the Appendix. The 
total covariance matrix is a sum of two terms: a sample 
covariance part accounting for the fact that different sea- 
sons of observation at different frequencies are observing 
the same CMB modes and that some of our cross spectra 
have common noise, and another part coming from the 
covariance of uncertainties in the determination of the 
beam profile. The covariance matrix is computed an- 
alytically and checked against Monte Carlo simulations 
described in Section [6l This total covari ance matrix is 
used in defining the like lihood function (jDunklcv et al.l 
l2013t iSievers et al.l 120131 ) when determining cosmologi- 
cal parameters. 

Along with the covariance matrices, we also generate 
bandpower window functions which convert a theoretical 
power spectrum into a band power: C^^ = J2e BbtCf^ . 
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Fig. 8. — Top: The IRIS 100 /jm map on the ACT-E strip (arbitrary color scale). Bottom: The equatorial dust mask based on the IRIS 
map flux cut as described in the text. The small box near right ascension of 00*" 14™ shows the "seagull" -like structure that is additionally 
masked out even though some pixels fall below the flux cut. 



Due to different geometry and noise properties of the 
ACT-E and ACT-S patches, two separate window func- 
tions Bm;Act-s and Bu-act-c are evahiated for the 
south and the equator. 

4.6. Combining Multi-Season Spectra 

For the purpose of parameter estimation, we keep the 
ACT-E and ACT-S spectra an d covariances for eac h sea- 
son separate in the hkchhood (jDunklcv et al.|[2013| ). For 
display purposes and for visual comparison with other 
datasets we combine the spectra from different seasons 
(separately for equator and south) using inverse variance 
weighting: 

jf^-l^{aAxPB);(aAxPB)^{aAx0B) 

f^AxB ^ l^a,/3y'^ )bb (O] 

b Y- /rs.^i-.{aAxpB):XaAxl3B) ' 

l^a,py^ )bb 

As discussed in iDunklcv ct al.l (j2013f ) the amplitudes 
of the Galactic cirrus contributions to the ACT-E and 
ACT-S maps are different. Therefore, before combin- 
ing the ACT-E and ACT-S spectra obtained above, we 
subtract the best-fit cirrus component (see 15.11 for more 
details) from the ACT-E and ACT-S spectra, and then 
combine them using inverse variance weighting. The 
multiple levels of cross-correlation used in computing the 
power spectrum help ensure that potential peculiarities 
in the observation that are located in time or space do 
not propagate to the final power spectrum. 

5. FOREGROUNDS 

In the 148 GHz and 218 GHz bands the main fore- 
grounds are emission from point sources and diffuse 
Galactic dust, which we treat with the application of 
masks as desc r ibed below . In the companio n papers 
IDunklcv et"all (j2013D and iSicvers ct all (|2013D we also 
consider and constrain contributions from thermal and 
kinetic Sunyaev Zaldovich effects, clustered and Poisson- 
like infrared point sources, radio sources, and a residual 
Galactic cirrus component. 

5.1. Galactic Dust 

We detect a significant contribution from Galactic cir- 
rus in our ACT-E maps, especially at 218 GHz. We em- 
ploy a two-step approach for dealing with Galactic cir- 
rus in the ACT-E maps using 100 jim du st maps from 
IRIS (jMiville-Deschenes fc Lagachd [20051) as the refer- 
ence. The first step is motived by the observation that 



most of the dust contamination in the equatorial power 
spectrum comes from the regions corresponding to bright 
clustered structures in the dust map. Therefore, we gen- 
erate a dust mask by identifying and setting to zero all 
pixels above a flux density of 5.44 MJy/sr as well as 
pixels that fall inside significantly clustered structures 
such as the "seagull" -like structure near right ascension 
of 00*^14™ shown inside the box in Fig. HI This dust mask 
is multiplied by the point source mask described below 
to generate the final mask that is applied to the data. 

The second step of the dust treatment is generating a 
model of the residual dust contamination after the ap- 
plication of this mask, and then to inform the parameter 
estimation pipeline with reasonable priors on this model 
(the residual model amplitude is fitted and marginal- 
ized over wh en constraining cosmo logical parameters, as 
described in IDunklcv et al.l I2013f). The mo d el is con- 
structed as follows. Following iHaiian et al.l (|2012f ) we 
perform a multicomponent fit to the auto-power spec- 
trum of the IRIS map after application of the dust mask 
described above. The components include a power law 
term for the residual Galactic cirrus, a Poisson shot noise 
term, a term representing the clustered component of in- 
frared emission, and a white noise term to describe the in- 
strument. The Galactic cirrus component is modeled as 
(^cirrus _ ^^j^^^g where the value of the power law 

index appears to be a good fit to Herschel observations o f 
cirrus ()Mivillc-Dcschcnes et al.l2010HBracco et al.l201lD . 
as well as the cross c orrelation between ACT and BLAST 
(jHaiian et al.ll2012D . and that between ACT and IRIS 
maps. This fitting procedure provides us with an esti- 
mate of the amplitude Acirrus separately for the ACT-E 
and the ACT-S map footprints. Next, we cross-correlate 
the ACT maps with the IRIS template to evaluate the 
dust coefficient A, = c^ct x iris /^misx mis ^^^^ 
quency and each sky region. Finally, the cirrus contri- 
bution to the ACT power spectrum can be expressed as 

^' = Acirrus^^ or cxprcsscd in terms of a rescaled 
amplitude at 4 = 3000: Be = e^Ci;/{2n) = Og (^/4)"° ^ 

where we have defined Og = ^chrus ^cr°'''/(27r). 

The various model parameters obtained from the fit- 
ting method above are displayed in Table O There is 
roughly twice as much dust in the ACT-E region as in 
ACT-S, but at e = 3000 and for 148 GHz it is less than 
3% of the CMB signal. These values represent a fre- 
quency scaling consistent w ith the early release results 
from the Planck satellite (jPlanck Collaboration et al] 



ACT Three-season Power Spectrum 



9 



TABLE 2 

Coefficients for Galactic cirrus model 



Region 


-^cirrus 

MJy2 


Frequency 
GHz 


MK/MJy 


ag 




ACT-S 


9.95 


148 
218 


8.65 
30.0 


0.4 
5.2 


ACT-E 


17.9 


148 
218 


8.65 
30.0 


0.8 
9.4 



^ we use a common dust coefficient for equator 
and south 



1201 It ), and can be compactly written in flux units as 
= agii/eo)-^-^!^,!^]/^^)^ with /3 = 1.8, i^' and 
the two frequencies being crossed, and i^o =148 GHz. 
Based on the scatter observed in these central values 
as well their variation depending on whether the clus- 
tering term is included in the fit, we adopt priors of 
ag, = 0.4±0.2 and age = 0.8±0.2 for ACT-S and ACT-E 
respectively. The fitted models and priors are illustrated 
in Fig. M 




Fig. 9. — Galactic cirrus contributions to the power spectrum 
modeled as a power law as described in the text, for each fre- 
quency and region of the sky. The curves correspond to the best 
fit amplitudes obtained by cross-correlating ACT maps with cirrus 
component of the IRIS maps, and the band around them repre- 
sent the un certainty adopted as a prior in the ACT likelihood as 
discussed in lDunklev et ahl I I2013I) . 



5.2. Point Sources 

Point sources are the main astrophysieal foreground for 
the 148 and 218 GHz bands. At these frequencies, point 
sources are typically either radio-loud AGN or dusty star- 
forming galaxies. Most of the bright sources are AGN, 
while most of the dusty star-forming galaxies lie below 
the detection threshold of our survey. Point sources must 
be identified and masked before the power spectrum is 
computed so as not to add power to the cosmological 



signal. We have id entified sources using a matched fi l- 
ter algorithm (e.g.. iTegmark fc de 01iveira-Costalll998D . 
We mask data within a 5' radius around all sources de- 
tected down to 15 mjy in either band. The residual 
power contributed to the power spectrum from unmasked 
sources below our detection threshold is expected to be 
2.9 ±0.4 jxK^ at I = 3000 (Gralla et al, in preparation). 
For details about the point source detection algorithm 
we used, and catalogs for the south 2008 148 GHz data, 
see [Marriage et al.l (j201lD . Catalogs for the remaining 
data set will be published in subsequent papers. 

6. SIMULATIONS 

We ran a set of Monte Carlo (MC) simulations in order 
to validate the analytic prescription for the uncertainties 
on the cross-season cross-frequency power spectra, and 
to investigate bandpower covariance and possible biases 
in the pipeline. As our map-making procedure is iter- 
ative, it is prohibitively expensive to run a large set of 
end-to-end simulations that would capture all aspects of 
the map-making pipeline, and the noise characteristics 
and corr elations in the actual data set. Instead, follow- 
ing ID 111 , we generate signal maps as Gaussian random 
realizations from a power spectrum, and add to each of 
them a realization of a Poisson point source population, 
and a simulated noise map generated from the observed 
noise-per-pixel in the data maps. The details of the im- 
plementation are essentially the same as in Section 4 of 
ID 111 with special care taken so that signal realizations 
are properly correlated across different seasons and foot- 
prints. For each season and each frequency, we generate 
960 signal-|- noise maps (four splits for each frequency), 
and for each realization we compute the power spectra in 
exactly the same way as we do for the data maps. From 
the large set of cross-power spectra obtained in this way 
we estimate the season-season covariance as well as the 
correlation between band powers. We find that in all 
cases, correlation between adjacent bins are insignificant 
at the 10% level. 

We evaluate the uncertainties in the band powers us- 
ing an analytic prescription described in Appendix lAl 
We verify the accuracy of these expressions by comparing 
the predicted error bars with the scatter of MC realiza- 
tions. For isotropic white noise realizations with uniform 
weights, our expressions for multi-season multi- frequency 
error bars are good to better than a percent. 

7. TEMPERATURE POWER SPECTRUM RESULTS 

Power spectra are computed following the procedure 
outlined in Section U separately for each region (south 
and equator) and for each season pair. The entire set of 
spectra along with their covariance is passed on to the 
likelihood code that forms the basis of parameter con- 
straints. Although combined spectra are not used in the 
actual analysis, in this section we discuss various combi- 
nations of power spectra for purposes of comparison and 
systematic tests. Note that ACT-S and ACT-E spectra 
cannot be trivially combined as residual Galactic cirrus 
contribution to the two regions are different. Therefore, 
we subtr act the best fit resi dual cirrus model (as dis- 
cussed in lDunklev et al.ll2013D from the estimated power 
spectra before combining ACT-E and ACT-S spectra. 
To simplify the presentation, all figures in this section 
portray dust-subtracted spectra. Another complication 
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TABLE 3 

Single frequency combined bandpowers provided for plotting purposes 

ONLY. 

% = 4(4 + l)Ci,/27r (hK2) 



148 GHz 148 GHz x 218 GHz 218 GHz 

I range central 1^ Bt Bt <^(Bh) ""(Sfc) 



540 - 


640 


590 


2267.4 


114.3 


640 - 


740 


690 


1760.2 


79.4 


740 - 


840 


790 


2411.2 


97.0 


840 - 


940 


890 


1962.4 


75.5 


940 - 


1040 


990 


1152.2 


42.8 


1040 - 


1140 


1090 


1208.7 


43.2 


1140 - 


1240 


1190 


1057.7 


36.0 


1240 - 


1340 


1290 


743.1 


25.7 


1340 - 


1440 


1390 


833.3 


27.3 


1440 - 


1540 


1490 


683.0 


22.0 


1540 - 


1640 


1590 


484.7 


16.4 


1640 - 


1740 


1690 


403.1 


13.3 


1740 - 


1840 


1790 


377.7 


12.3 


1840 - 


1940 


1890 


266.8 


9.3 


1940 - 


2040 


1990 


236.5 


8.7 


2040 - 


2140 


2090 


229.2 


8.1 


2140 - 


2340 


2240 


150.2 


4.3 


2340 - 


2540 


2440 


109.2 


3.5 


2540 - 


2740 


2640 


75.0 


3.0 


2740 - 


2940 


2840 


63.7 


2.9 


2940 - 


3340 


3140 


43.9 


1.9 


3340 - 


3740 


3540 


35.7 


2.0 


3740 - 


4140 


3940 


36.1 


2.4 


4140 - 


4540 


4340 


34.4 


2.8 


4540 - 


4940 


4740 


37.3 


3.4 


4940 - 


5840 


5390 


50.5 


3.2 


5840 - 


6740 


6290 


59.5 


5.2 


6740 - 


7640 


7190 


81.6 


8.7 


7640 - 


8540 


8090 


131.4 


14.6 


8540 - 


9440 


8990 


133.0 


25.5 



arises due to the different geometries and masking pat- 
tern of the ACT-E and ACT-S maps, which eause the 
theoretical bandpowers for these regions to be in prin- 
ciple different, although the actual differences are small. 
Also, due to subtle variations of the beam profile from 
one season to another, the beam uncertainties in individ- 
ual season spectra are slightly different. All these sub- 
tleties prompted the separate treatment of power spectra 
in the likelihood. In this section, we neglect these sub- 
tleties and combine spectra, with inverse variance weight, 
across season pairs and regions of the sky. We warn the 
reader hat such combinations are for visualization pur- 
poses only. Fig. [TUl shows the ACT-E and ACT-S spectra 
combined across the different observing seasons, along 
with their corresponding theoretical band powers. The 
power spectrum combined across all seasons and across 
the ACT-E and ACT-S strips is displayed in Fig. El 
The corresponding band power values and uncertainties 
are tabulated in Table [31 These plots portray how our 
pipeline is able to produce an estimate of the power spec- 
trum over the entire multipole range of 500—10000. Over 
the multipole range of ^ ~ 500—2500 these spectra clearly 
show the Silk damping tail of the CMB power spectrum, 
while on smaller angular scales {i ^ 2500 — 10000) a clear 
excess from the frequency-dependent Sunyaev Zcl'dovich 
effects and extragalactic foregrounds (radio and infrared 
point sources) is cl early visible (these c ontributions are 
further discussed in lDunklev et al.ll2013f ). 

Finally, we display, in Fig. [121 the state of the art 
in CMB temperature power spectrum measurements 
down to the damping tail where we plot the WMAP 











494.0 


15.7 


551.2 


30.0 


400.6 


12.6 


458.5 


26.9 


369.7 


11.5 


408.3 


23.9 


272.7 


9.3 


327.0 


20.3 


261.3 


8.8 


320.2 


20.5 


226.6 


7.7 


274.7 


17.9 


168.6 


4.6 


238.7 


12.3 


128.1 


3.8 


199.8 


10.7 


97.5 


3.3 


181.6 


9.8 


86.0 


3.2 


181.9 


9.6 


69.4 


2.1 


167.3 


6.4 


65.2 


2.1 


183.6 


6.5 


68.3 


2.3 


211.4 


7.3 


75.5 


2.7 


245.7 


8.3 


93.8 


3.3 


286.8 


9.8 


109.8 


3.0 


355.8 


9.5 


149.6 


4.6 


478.0 


14.6 


177.9 


6.6 


564.4 


19.9 


240.1 


10.7 


753.1 


29.6 


265.3 


16.6 


878.7 


44.3 



7-year spectrum, the inverse variance combined ACT- 
E+A CT-S spectrum, a nd the recent SPT power spec- 
trum (|Storv et al.ll2012D . 

7.1. Power spectrum with alternative binning 

As discussed in Section [431 the choice of large bins for 
our main power spectrum result was motivated by the 
need for keeping the bandpowers minimally correlated. 
It is of interest, however, to ask how the spectrum would 
have looked with smaller bi ns of w idth = 50 over the 
damping tail, as was done in ID 111 Such a result is shown 
in Fig. [T3l Note that the first through the eighth peak 
of the CMB can be clearly seen with this binning. We 
do not pursue this binning any further for the aforemen- 
tioned reasons. 

7.2. Derived CMB-only power spectrum 

The ACT-E and ACT-S spectra shown in Fig. \M in- 
clude the primary CMB signal as well as power from fore- 
grounds and SZ. We show the estimated primary CMB 
spectr um from ACT in Fig.[T31 derived in lDunklev et al.l 
()2013[ ). There, the multi- frequency spectra are used to 
estimate the CMB in bandpowers for ACT-E and ACT-S, 
simultaneously with the SZ and foregrounds components. 
The CMB spectra for ACT-E and ACT-S are then coad- 
ded for display. No assumptions are made about the 
cosmological model, only that the CMB is blackbody. 
Using the multi-frequency data to separate components, 
the CMB power can be recovered out to mrdtipoles of 
i ~ 3500. 
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Fig. 10. — Combined multi-season power spectra for the ACT-E Strip (upper panels) a nd the ACT-S Strip (lower panels). For the ACT-E 
and ACT-S data points, the corresponding best fit residual cirrus model (as discussed in lDunklev et ani2013l) has been subtracted. The left 
hand panel shows a linear scale zoomed-in version of the spectrum with an scaling to emphasize the higher order acoustic features. The 
lines s how the binned versio n of the best fit model for each frequency pair including CMB secondaries and foregrounds from our companion 
paper IDunklev et al.l l |2013l ) . The right panel shows the entire range of the com puted spectrum on a l og-linear scale with the conventional 
1(1 -\- 1) scaling. The lines show the unbinned version of the best fit model from IDunklev et al.l II2013I ). 



7.3. Systematic Tests 

In order to check for systematics in the map-making 
and power spectrum estimation pipelines, we perform 
various tests on the data. These are constructed such 
that the sky signal cancels between the various splits of 
the data, and only systematic effects remain. We test 
that the power spectrum obtained is the same in each 
season, in all time splits, from different parts of the ar- 
ray, with and without data near the telescope turnaround 
points, from different directions in Fourier space, and for 
different regions of the sky. The statistics from a subset 
of these tests are summarized in Tables lU [5l [6l and [T] 
Given that the spectra are computed individually and 
then included in the likelihood with the full covariance 
of the different frequencies and seasons, we compute the 
null tests on each subset of data, both for ACT-E and 
ACT-S, and for different seasons. 

7.3.1. Cross season nulls 



TABLE 4 

Null test values for the season consistency 

TESTS PERFORMED ON THE ACT DATA. THE PROBABILITY 
TO EXCEED (PTE) THE IS SHOWN IN PARENTHESES. 



Frequency 


Region 


Seasons Seasons Seasons 


dof 






2008-2009 2009-2010 2008-2010 




148 GHz 


South 


32.2 (0.36) 30.7 (0.43) 35.0 (0.24) 


30 




Equator 


39.7 (0.11) 


30 


220 GHz 


South 


27.7 (0.11) 15.9 (0.72) 21.5 (0.37) 


20 




Equator 


24.2 (0.23) 


20 



First, we test the year-to-year consistency of power 
spectra. In order to account for differences in the ACT 
beam from one observing season to the next, wc con- 
volve the map from one season with the beam profile of 
the other season being differenced, so that each map ef- 
fectively has the same beam transfer function. Then we 
difference the corresponding splits from seasons si and 

Ar(n)^[r;(n)-r;(n)], (3) 

where i = 1.2,3,4 represent the split index. The pixel 
weight map W corresponding to these difference splits is 
computed as: 

W-'^W-' + W~\ (4) 

where Wg-^ is the total iVobs map for season si etc. As 
with the other null tests, the azimuthal weighting is com- 
puted using the weights from the full data spectrum run. 
Figure [T3] shows these test for the ACT data, while the 
individual values for the tests are summarized in Ta- 
ble m We find all spectra computed in this way to be 
consistent with null. 

7.3.2. Split Nulls 

As discussed in Section [H the data in each season are 
separated into four splits in such a way that the detector 
noise is independent from one split to another. There- 
fore, the difference between any two splits should be con- 
sistent with noise and the signal should subtract away. 
We test this by generating difference maps from each 
pair, and computing the two-way cross spectra from in- 
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2000 4000 6000 8000 10000 

Multipole I 

Fig. 11. — Multifrequency power spectra combined across all seasons and the ACT-E a nd ACT-S regions. Be fore combining the ACT-E 
and ACT-S data points, the corresponding best fit residual cirrus model (as discussed in lDunklev et ani2013 ) has been subtracted. The 
upper panel shows the £ = 500 — 4000 portion of the power spectrum on a linear scale with an £^ scaling to emphasize the higher order 
acoustic peaks. The lower panel shows the entire range of the computed spectra with the l{t. + 1) scaling. The lines in either case sho w the 
best fit models for each frequency pair including CMB secondaries and foregrounds from our companion pap er IDunklev et al.l l|2013l ^ . The 
grey data points represent the power spectrum from the WMAP seven-year data release l ILarson ct al.|[2"011l) . 
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Fig. 12. — State of the art of CMB temperature power spectrum m easurements from the WMAP 9-year data release l IBennctt et al.ir2012l : 
IHinshaw cFaIll2012l) . the South Pole Telescope llStorv ct al.ll2012l) and AC T (this work). The solid line shows the best fit model to the 
ACT 148 GHz data combined with WMAP 7-vear data (ILarsorrct al.ll2011l ). The dashed line shows the CMB-only component of the same 
best fit model. Although we compute the power spectrum down to £ = 200, we do not use data below £ = 540 in the analysis. 



TABLE 5 

^2 (^p'pg-j VALUES FOR THE TOD SPLIT NULL TESTS PERFORMED ON THE ACT DATA. 



Frequency 


Region 


Season 


TOD 


dof 








(l-2)x(3-4) (l-3)x(2-4) (l-4)x(2-3) 




148 GHz 


South 


2008 
2009 
2010 


19.7 (0.92) 37.7 (0.16) 37.1 (0.17) 
30.3 (0.45) 31.8 (0.38) 22.6 (0.83) 
35.7 (0.22) 28.6 (0.54) 21.5 (0.87) 


30 
30 
30 




Equator 


2009 
2010 


33.9 (0.29) 26.5 (0.65) 40.9 (0.09) 
34.6 (0.26) 35.6 (0.22) 24.0 (0.77) 


30 
30 


220 GHz 


South 


2008 
2009 
2010 


33.1 (0.03) 28.2 (0.10) 15.3 (0.76) 
14.4 (0.81) 11.3 (0.94) 14.8 (0.79) 
8.8 (0.99) 16.0 (0.72) 21.0 (0.40) 


20 
20 
20 




Equator 


2009 
2010 


24.9 (0.21) 19.3 (0.50) 13.3 (0.87) 
11.8 (0.92) 16.3 (0.70) 14.0 (0.83) 


20 
20 
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Fig. 13. — Combined ACT-E + ACT-S 148 G Hz power spectrum computed with alternate binning shown alongside the WMAP 9-year 
data IIBennett et al.|[2012l : [Hinshaw et al.|[20T^ ). Note that with these smaller bins, the contours of the first seven acoustic peaks of the 
CMB power spectrum can be clearly seen. The bandpowers are significantly correlated at this bin size, and a precise estimate of the bin 
to bin correlation is computationally pro hibitively costly. The solid line shows the best fit model to the ACT 148 GHz data combined with 
WMAP 7- year data IILarson et al.ll2011l) . 



dependent pairs of difference maps, e.g.: 

ri2(n) = [T^{h) -T^h)]/2 

T'\h)^[T'ih)~T\n)]/2. (5) 

The difference maps are expected to contain noise but no 
residual signal. We estimate the cross-spectrum of the 
difference maps, Cb = (T^'^T^'^^ , and the other two per- 
mutations of the differences ((T^T^^) and {T^'^T'^'^)) . 
These difference maps are downweignted by the same 
weight maps used to construct the full power spectrum. 
Similarly, the azimuthal weights are borrowed from the 
full data spectrum run. The three difference spectra are 
shown in FigureHHfor the 148 GHz ACT-S data set. The 
statistics corresponding to this test are shown in Table [5l 
The spectra are found to be consistent with a null signal, 
as expected. 

7.3.3. In/Out nulls 

In order to test for systematic detector asymmetries, 
we make a map using data from detectors from the inner 
region of the array, and another map from detectors along 
the edges, and compute the differences between the two 
maps: 



Tf{n)]/2 
T/(n)]/2 



(6) 



where the i and o label the inner and outer parts of the 
detector array respectively. The full set of values are 
summarized in Table [51 In general we see no trend for 
differences as a function of detector position; the null 
tests arc consistent with no signal. 

7.3.4. Turnarounds 



Another null test, based on cutting out data around 
telescope turnaround probes the consistency of data 
taken with the telescope accelerating as it reverses direc- 
tion at the ends of the scan (turnarounds). In the maps 
used for the power spectrum estimation, the data taken 
during the turnaround is included. We test for any arti- 
facts generated by the acceleration at turnaround by tak- 
ing the difference of maps with and without turnaround 
data. Maps are made cutting data near the turnarounds, 
amounting to removing w 10% of the total data. This 
loss of data affects the two sky regions differently. In the 
southern patches, the loss of data is uniform and leads to 
a slight increase in striping in the maps, whereas in the 
equatorial patches, removing the turnarounds removes 
data at the upper and lower edges of the maps. Hence, 
for these tests we compute the differences using an equa- 
torial region which is slightly narrower in the declination 
direction (2? 7 as opposed to 3dcg wide). Two difference 
maps are made by pairing one split of the standard map 
with a different split of the new maps (we avoid differ- 
encing the same splits as they have very similar noise 
structure), and a two-way cross-power spectrum is pro- 
duced. Any artifact due to the turnaround would be left 
in these difference maps and might produce excess power. 
We compute the turnaround cuts as a function of season, 
frequency range and area on the sky. The reduced val- 
ues are summarized in Table [7l Again, we find that the 
difference maps have spectra consistent with no signal. 

7.3.5. Distribution 

While the null tests are performed for different subsets 
of the data, we combine the statistics from the null tests 
together to test for consistency globally. We restrict the 
range of the 218 GHz spectrum to be ^ > 1500, hence the 
218 GHz spectrum contains 20 degrees of freedom, while 
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Fig. 14. — The CM B power spectrum estimated from ACT, shown with the spectrum from the WMAP 9-year data l IBennett et al.ir2012l ; 
IHinshaw tiTSiXWvh . The errors include uncertainty due to foreground and SZ emissi on, as well as the rela tive calibration of the 148 
and 218 GHz channels, and beam uncertainty. The full covariance matrix is derived in lDunklev et aTl II2013I ). The solid line shows the 
CMB-only component of the best fit model for the ACT data combined with the WMAP 7-year data. 



TABLE 6 

Null test (PTE) values for the inner vs outer detectors. 



Frequency 


Region 


Season 


In/Out 


dof 








(lo-2i)x(3o-4i) 


(lo-3i)x(2o-4i) 


(lo-4i)x(2o-3i) 




148 GHz 


South 


2008 


26.2 (0.66) 


28.2 (0.56) 


26.2 (0.67) 


30 






2009 


37.7 (0.16) 


17.3 (0.97) 


27.3 (0.61) 


30 






2010 


34.7 (0.26) 


27.3 (0.61) 


25.6 (0.70) 


30 




Equator 


2009 


32.0 (0.36) 


24.8 (0.74) 


26.2 (0.66) 


30 






2010 


32.2 (0.36) 


36.3 (0.20) 


38.6 (0.13) 


30 


218 GHz 


South 


2008 


25.3 (0.19) 


23.9 (0.24) 


21.0 (0.40) 


20 






2009 


13.7 (0.85) 


12.3 (0.91) 


25.8 (0.17) 


20 






2010 


13.8 (0.84) 


15.0 (0.78) 


26.3 (0.16) 


20 




Equator 


2009 


23.4 (0.27) 


22.6 (0.31) 


27.9 (0.11) 


20 






2010 


9.9 (0.97) 


25.3 (0.19) 


9.2 (0.98) 


20 
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Fig. 15. — Cross season null tost for the ACT data. The top row 
illustrate the 148 GHz cross season null tests for ACT-E (top left) 
and ACT-S (top right), while the bottom row show the 218 GHz 
eross seaso n nu lls. The values for the fit are presented in Table|4] 
and Figure [T71 
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Fig. 17. — The reduced x values for all null tests. The blue 
histogram is computed for the values from the 218 GHz null 
tests, while the purple histogram shows the same null tests for 
the 148 GHz maps. The black dashed and dot-dashed lines show 
the theoretical distributions for 20 (AR2) and 30 (ARl) degrees 
of freedom respectively, normalized to match the frequency of the 
histograms. The x^ values presented here are given in Tables |4l (5] 
[6] and [3 
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Fig. 16.— TOD null test for the 148 GHz Southern strip, from 
2008 (top panel) to 2010 (bottom panel). For each year, three 
TOD nulls are created from the combinations described in Eq. [5] 
The x^ values for the null test are summarised in Table [5| 



the 148 GHz spectrum containts 30 degrees of freedom. 
We show the distribution of values and the theoreti- 
cal distribution for the two cases in Figure [T71 This 
shows that the null tests are broadly consistent with be- 
ing drawn from a distribution for the number of de- 



7.3.6. Isotropy 



e=\ogM£ + lf Cj2Tr /{I iiK^ ) I) 




1/ 1000 



Fig. 18. — The two dimensional 148 GHz cross-power spectrum 
co-added across the ACT-E patches and seasons. For £ < 2500 
smoothing by a small kernel has been applied. The acoustic fea- 
tures in the power spectrum are clearly visible. At £ > 2500, where 
the instrument noise dominates, we display the raw spectrum. The 
black lines represent the CMB-only theory and have been plotted 
to guide the eye. 



We test the isotropy of the power spectrum by estimat- 
ing the power as a function of phase 9 = arctan(€j,/f^). 
We compute the inverse-noise-weighted two-dimensional 
spectrum co-added across patches and seasons for the 
ACT-E region. We show the mean two-dimensional 
cross-power pseudo spectrum in Figure [THl The spec- 
trum is symmetric for £ to —£, as it is for any real valued 
map. To quantify any anisotropy, the power averaged 
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TABLE 7 

Null test values for turnarounds 



Fr c Quc 11 c y 


Region 


Season 


Turnarounds 


dof 








(lt-2nt)x(3t-4nt) 


(lt-3nt)x(2t-4nt) 


(lt-4nt)x(2t-3nt) 




148 GHz 


South 


2008 


21.9 (0.85) 


24.8 (0.74) 


34.6 (0.26) 


30 






2009 


28.3 (0.56) 


35.3 (0.23) 


26.7 (0.64) 


30 






2010 


30.5 (0.44) 


30.4 (0.44) 


35.0 (0.24) 


30 




Equator 


2009 


27.4 (0.60) 


29.3 (0.50) 


34.0 (0.28) 


30 






2010 


35.1 (0.24) 


43.9 (0.05) 


21.9 (0.86) 


30 


218 GHz 


South 


2008 


25.9 (0.17) 


25.3 (0.19) 


20.6 (0.42) 


20 






2009 


15.8 (0.73) 


13.0 (0.88) 


14.3 (0.82) 


20 






2010 


11.2 (0.94) 


11.6 (0.93) 


21.1 (0.4) 


20 




Equator 


2009 


31.3 (0.06) 


19.4 (0.5) 


16.6 (0.68) 


20 






2010 


14.1 (0.82) 


24.3 (023) 


16.7 (0.67) 


20 



over all multipoles in the range 200 < £ < 10000 is com- 
puted in wedges of 6' = 20° , and compared to the mean of 
the entire annulus. No significant deviation from isotropy 
is detected using this method. We find that this result 
holds for ACT-S and 218 GHz maps. 

7.4. Consistency of ACT-E, ACT-S, and SPT spectra 

As mentioned above, due to the difference in geometry 
of the equatorial vs. southern patches, the band power 
binning functions for ACT-E and ACT-S are slightly dif- 
ferent leading small differences in t he binned version o f 
the best fit ACT + WMAP7 model (jSievers et al.|[20Tl . 
Therefore to test the consistency between ACT-E and 
ACT-S spectra we check for the nullity of the residuals 
from their corresponding binned best fit model. We also 
consider the consis tency of the SPT band powers from 
iKeisler et al.l ()2011[ ). Care must be taken while comput- 
ing the residual for the SPT spectrum as the point source 
masking threshold corresponding to that spectrum was 
different from that of ACT. To correct for this, we adjust 
the Poisson point source componen t of our best fi t mode l 
to match the masking level used in IKeisler et al.l ()2011l ). 
The results are shown in Fig. [191 clearly indicating that 
the three spectra are consistent with null, and therefore 
with each other. 

The suite of consistency tests performed here show that 
our reported spectrum passes a wide range of checks for 
systematic errors in time, detector-space, map-space, and 
^-space. 

8. GRAVITATIONAL LENSING ANALYSIS 

Large-scale structure gravitationally deflects the CMB 
radiation as it passes through the universe, thereby defin- 
ing a lensing deflection field d that remaps the observed 
CMB temperature sky: T(n) = Tunicnscd(n-l-d). Lensing 
distorts the small-scale CMB anisotropics, thus modify- 
ing their statistical properties, with the lensing deflection 
locally breaking statistical isotropy and correlating for- 
merly independent temperature Fourier modes (or more 
intuitively, correlating the CMB temperature with its 
gradient). The lensed small-scale CMB thus contains not 
only information about the universe at the last-scattering 
surface (z ~ 1100), but also encodes information about 
the cosmic mass distribution at later times and lower red- 
shifts. Using an estimator quadratic in temperature that 
measures the lensing-induced change in the statistics of 
the CMB fluctuations from the correlation of formerly 
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Fig. 19. — The residual power speetra obtained after subtract- 
ing the best fit ACT -I- WMAP model from the ACT-E, ACT-S, 
and SPT power spectrum (IKeisler et al.ll2011l ). For the SPT resid- 
ual the Poisson point source component of the best fit model is 
adjusted to reflect the difference in point source masking levels 
between ACT and SPT. The residuals are all null showing the con- 
sistency of these spectra. 

independent Fourier modes, one can construct a noisy 
estimate of the CMB lensing convergence, k = — |V • d, 
which is a measure of the projected matter density out to 
high redshifts. The power spectrum of the CMB lensing 
convergence can be simply obtained from this estimated 
lensing convergence field, though biases arising from in- 
strument and cosmic-variance-induced noise as well as 
from higher order corrections to the estimator must be 
subtracted. Equivalently the estimation of the lensing 
power spectrum can be regarded as a measurement of 
a lensing-induced non-Gaussianity in the CMB temper- 
ature four-point correlation function. Measurements of 
the lensing power spectrum can place strong constraints 
on the properties of dark energy and neutrinos, and also 
serve as a valuable test of the ACDM prediction for struc- 
ture growth and geometry at redshifts z ~ 1 — 3. 

CMB lensing science has made great advances in re- 
cent years. Signatures of CMB lensing were first ob- 
served in correlations of the CMB with large-scale struc- 
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Fig. 20. — CMB convergence power spectrum reconstructed from 
the ACT-E strip temperature data. The enhanced effective depth 
of the coadded ACT-E map (~ 18 /jk-ar cmin) compared t o its 
previous version (~ 23 /^k-arcmin) used in lDas et al.l l|2011al ) led 
to an improved detection significance. 



ture (jSmith ct al.ll2007l : iHiTata et al.ll2008D . The lensing 
power spectrum was fir st detected at Aa significance by 
the ACT collaboration pas et al.ll2011aD . The first evi- 
dence for dark energy from the CMB alone was also ob- 
tained from this measurement bv lSherwin et al.ll2011l A 
more significant detection (at 6.3 a) of the lensing power 
spectrum and more precise cosmological constraints were 
reported by the SPT co llaboration the following year 
(jvan Engelen et al.|[20T^ . 

In this work we present a measurement of the lensing 
power spectrum from the improved ACT maps on the 
ACT-E strip. Our new measuremei it of lensing essen - 
tially use s the s ame methodology as iDas et al.l (j2011af ). 
hereafter IDS 111 . Lensing is measured using a quadratic 
estimator in temperature; the power spectrum of the 
CMB lensing convergence is thus a temperature four- 
point functi on wit h the filtering and normalization ob- 
tained as in Ids 111 . The bias is calculated in three steps 
to make the calculation as robust and model-independent 
as possible. First, we simulate the bias that is present in 
the absence of any mode-correlations, removing existing 
correlations by randomizing the phase s of each Fourier 
mode of the temperature field (jPSllf ). This is mathe- 
matically equivalent to estimating the so-called N{Q) bias 
from the pscudo-Cf power spectrum measured from this 
map. The use of measurements rather than simulations 
makes the bias calculation more robust. Next, we simu- 
late additional small biases from finite-map effects, and 
anisotropic and inhomogeneous noise by performing lens- 
ing power spectrum estimation on unlensed simulated 
CMB maps with realistic noise characteristics. After ac- 
counting for these two sources of bias, the input lensing 
power spectrum is recovered fairly accurately in simula- 
tions, but there is a small difference between recovered 
and input power spectra due to higher order corrections. 
We adopt the small, simulated difference between the 
recovered and true power spectra as an additional bias 
and subtract it from the biased lensing power spectrum 
estimated from data. 

Systematic contamination of the estimato r by th e SZ 
signal, IR and radio s ources was est i mated in lDSllI using 
the simulations from iSehgal ct afl ()2010D . The authors 
found that, with the ACT lensing pipeline as used in this 



work, the contamination is smaller than the signal by two 
orders of magnitude and can thus be neglected. This re- 
sult appears well-motivated for two reasons, which also 
apply to our current analysis: first, in the analysis we 
only use the signal-dominated scales below £ = 2300, at 
which SZ. IR and radio power are subdominant; second, 
by using the data to estimate the bias, our estimator au- 
tomatically subtracts the Gaussian part of the contami- 
nation, so that only a very small non-Gaussian residual 
remains. The previous ACT contamination estimates are 
not strictly applicable to the new lensing estimate, be- 
cause the filters used here contain somewhat lower noise, 
and thus admit slightly more signal at higher .^s; how- 
ever, estimates by the SPT collaboration (van Engelen 
et al. 2012) with similar noise levels and filters also find 
negligible contamination. The contamination levels in 
our analysis arc thus also expected to be negligible. 

The measured CMB lensing power spectrum, detected 
at 4.6cr, is shown in Fig. [201 along with a theory 
curve showing the convergence power spectrum for a 
fiducial ACDM model defined by the parameter set 
(f7fc,a„,f7A,/i,n„cr8) = (0.044, 0.264, 0.736, 0.71, 0.96, 
0.80). Constraining the conventional lensing parameter 
A]^ that rescales the fiducial convergence power spectrum 
{C^^ AlC^") we obtain Al = 1.06 ± 0.23. The data 
arc thus a good fit to the AC DM pr ediction for the ampli- 
tude of CMB lensing. As in IDS 111 , we find the spectrum 
to have Gaussian errors, uncorrelated between bins. For 
some parameter runs, the lensing power spectrum infor- 
mation is added to t he CMB power spectrum information 
(|Sievers et alllMl ). 

9. CONCLUSIONS 

We have derived the power spectrum of microwave sky 
maps at 148 GHz and 218 GHz produced by the At- 
acama Cosmology Telescope experiment, such as those 
displayed in Fig. 5. The power spectra cover a range 
of angular scales spanning nearly a factor of 20, rang- 
ing from around 0.35 degrees {£ = 590) to a little over 
one arcminutc {£ ~ 8900). The maps are high qual- 
ity, and in principle extracting the power spectrum is 
a simple matter. A host of practical considerations, 
along with the precision supported by the data, make 
estimation of the power spectrum challenging. This pa- 
per summarizes algorithms and techniques for handling 
the particular shapes of our maps, point source contam- 
ination, the steepness of the power spectrum, significant 
features due to galactic dust emission, spatially varying 
noise levels, and calibration. In addition to the resulting 
power spectra, we also display numerous null tests on 
the data. These tests, along with results from simulated 
maps, make a strong case that any systematic errors in 
our power spectra are below the level of statistical error. 

The ACT power spectra are consistent with those mea- 
sured by the South Pole Telescope collaboration, as are 
the underlying maps in a region of overlapping sky cov- 
erage. Given how small the signals are and how many 
sources of error must be tamed to measure them, consis- 
tent results represent a substantial experimental achieve- 
ment. 

The temperature power spectrum measurements dis- 
played in Figure 12 represent the culmination of a two- 
decade quest, since the first large- angle po wer mcasure- 
ments were made by the COBE satellite (|Smoot et al.l 
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[l99l) . It was soon realized that for inflationary cosmo- 
logical models, the substantial structure in the microwave 
background temperature angular power spectrum due to 
coherent acoustic oscillations in the early universe would 
allow pre cise constraints on th e basic properties of the 
universe (jJungman series of innovative 

and increasingly sensitive experiments then gradually un- 
veiled the power spectrum. With the definitive measure- 
men ts down to quarter-d egree scales by the WMAP satel- 
lite (|Bennett et al. 1120 121 ) and the precise arc minute-scale 
meas urements by ACT (this work) and SPT (jStorv et alJ 
l2012t ) along with the results anticipated from the Planck 
satellite, this particular route to cosmological insight is 
approaching a highly refined state. 

A new frontier in microwave background experi- 
ments will likely be detailed Icnsing m aps from high- 
resolution polarization meas urements (jNiemack et al.l 
l201d : lAustermann et aT]l2012| ). which have the prospect 
of co nstraining dark en ergy and modified gravity 
(e.g., iDas k Linde^ |2012[) . The lensing power spec- 
trum measurements presented here and by SPT 



(|van Engelen et ahl I2012D are the first steps along this 
new path. 

This work was supported by the U.S. National Sci- 
ence Foundation through awards AST-0408698 and AST- 
0965625 for the ACT project, as well as awards PHY- 
0855887 and PHY-1214379. Funding was also provided 
by Princeton University, the University of Pennsylva- 
nia, and a Canada Foundation for Innovation (CFI) 
award to UBC. ACT operates in the Parque Astronomico 
Atacama in northern Chile under the auspices of the 
Comision Nacional dc Investigacion Cientffica y Tec- 
nologica dc Chile (CONICYT). Computations were per- 
formed on the GPC supercomputer at the SciNet HPC 
Consortium. SciNet is funded by the CFI under the 
auspices of Compute Canada, the Government of On- 
tario, the Ontario Research Fund - Research Excellence; 
and the University of Toronto. SD acknowledges sup- 
port from the David Schramm Fellowship at Argonne 
National Laboratory and the Berkeley Center for Cos- 
mological Physics fellowship. RD acknowledges support 
from FONDECYT and BASAL grants. 



APPENDIX 



ANALYTIC ERRORBARS 



Here we derive an analytic expression for the expected error bars on the cross-frequency mulitseason cross-power 
spectrum. We denote the freq uenci es with uppercase A, B, C, D, the seasons with a, $, 7, r, and the sub-season data 
splits with i, j, k,l. Following ID 111 , the covariance the of cross-power spectrum is defined as: 



e 



-,(CxD). 



bb 



which expands as 



(Al) 



e 



{aAxPB)-XlCxTD) 
bb 



^ ^ \ " \ " (\ / rj^*iArj^3Brj^*kCrT^lD 



^iiAx]B) \ /^{kCxlD) 
^b,al3 I X^b.-yr 



X(l - SijSap){l - SmSjt)- 



(A2) 



The Kronecker deltas remove the auto power spectra, and any same-split, same-season cross- frequency spectra. The 
general normalization is 



i^j.k.l 



(A3) 



Applying Wick's Theorem, we have 
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Therefore, q(^^^^^)'(t'^^'^^) expands to 
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which after some algebra reduces to 
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Therefore, with A B C and a ^ (3 we have. 
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nd nd{nd - 1) 
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- 1) 



2Ct + 2- 
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Vb 



1 

Vb 

1 

Vb 
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CbN, 



aa.AA 



nd 



(A8) 

(A9) 

(AlO) 
(All) 
(A12) 

(A13) 
(A14) 
(A15) 



INCLUDING BEAM COVARIANCE IN THE ESTIMATE 

We can also include the effect of an uncertainty on the beam when combining data between different season. We 
consider the covariance of the power spectrum for two season pairs: i x js and k x I, with beam window functions 



,1X3 



[ij, — BlBl and Wf^,^ = B^iB^,, respectively. The measured windows function is given by ~ Wb + Swb- This 



error propagates to the power spectrum covariance as 



..ixjf^ixj kxl/^kxl r^ixj ^kxl 

I ' ixi_.fcxi\ b' / 



IX 



ixj „,,kxl 



The error on the window function is related to the error of the beam by 

Swi:'' =5BlBl+BlSBl. 



(Bl) 



(B2) 
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